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Mh. must be subtracted from the nuclear phase shifts obtained from the analysis of low 
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magnetic potential and an effective hadronic potential. We compare our results with 
those of previous calculations and estimate the uncertainties in the corrections. 
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1 Introduction 



The aim of this paper is to present a careful recalculation of the electromag- 
netic corrections which have to be applied in the analysis of low energy vr+p 
elastic scattering data in order to recover the hadronic phase shifts. 

For these electromagnetic effects in ttN scattering {n'^p elastic and 7r~p elastic 
and charge exchange) two different approaches were developed: the dispersion 
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theory method and a method using potentials for the hadronic and electro- 
magnetic interactions in a relativised Schrodinger equation (RSE for short). 
By this we mean an extension of the Schrodinger equation to the relativistic 
domain. 

The dispersion theory approach was initiated by Sauter [1,2]. It was then 
extended by the NORDITA group. The final results for their corrections are 
given in Tromborg et al. [3]. In Ref. [4] the same authors present a detailed 
exposition of the underlying ideas and point out their limitations. 

The potential theory approach was initiated by Oades and Rasche [5,6] and 
applied in detail by Zimmermann [7,8] to an analysis of the then existing data. 
His final results appear in the Landolt-Bornstein tables on tiN scattering [9]. 

Both methods have their shortcomings. In the dispersion theory approach 
the only contribution from t- and -u-channel exchange which is taken into 
account is that from -u-channel nucleon exchange. The rest of the unphysical 
contributions are omitted. Ref. [4] states: 'We are not able to determine or even 
estimate this term, which is of course a serious drawback of our method'. For 
this reason, Ref. [4] makes it clear that the results of an analysis of data for 
the three measured reactions which uses the corrections given there cannot 
be used to draw any conclusion about isospin symmetry violation. It seems 
that the fact that a dispersion theory method cannot give absolute values of 
the electromagnetic corrections has not been sufficiently recognised. This is 
an unsatisfactory situation, particularly for the analysis of the large amount 
of data at low energies. 

In the potential theory approach of Refs.[7,8] hadronic potentials were con- 
structed via RSEs to reproduce the hadronic phase shifts. These potentials 
then were used in RSEs, containing also the electromagnetic potential, to cal- 
culate the electromagnetic corrections. The idea is that by this addition of an 
effective short range hadronic potential and the long range electromagnetic 
potential we can estimate the corrections reliably. The weak point here is the 
use of a hadronic potential in a RSE. But it should be noted that 
(i) the hadronic potentials are used only to calculate the corrections, not to 
calculate the hadronic quantities themselves; 

(ii) the potential model is the only method that can give a well defined sepa- 
ration between hadronic quantities and their electromagnetic corrections. 

It was indeed possible (Ref. [8]) to analyse the best experimental data of that 
time in an isospin invariant way, using energy dependent hadronic potentials. 
The main reasons for recalculating the electromagnetic corrections in the po- 
tential model approach are the following. 

(i) The results need to be extended to lower energies, where there is now a 
large amount of new data obtained at pion factories. 



(ii) The slight effect which a change of the hadronic input has on the electro- 
magnetic corrections needs to be controlled. This has not been done by the 
NORDITA group. 

(iii) The effect of some fine details of the electromagnetic interaction, which 
might become important at the present precision of the experimental data, 
needs to be checked. This also is not contained in the NORDITA results, 
(iv) Instead of the energy dependent potentials used in Refs.[7,8] we want to 
use energy independent hadronic potentials to calculate the electromagnetic 
corrections for pion laboratory kinetic energy T^^ < 100 MeV. 

In this paper we treat only the single channel problem of 7r~^p elastic scattering. 
It enables us to describe the formalism carefully and to provide a survey of 
the relative importance of the contributions to the corrections. In the next 
paper we extend the formalism to the two-channel {7i~p, Ti'^n) case and give 
the results for the electromagnetic corrections which need to be applied in the 
analysis of low energy n^p elastic and charge exchange scattering data. 

In Section 2 we give a detailed description of our model for the electromagnetic 
corrections to 7r~^p scattering. In Section 3 we explain the method of evaluation 
of the corrections in this model and in Section 4 we give the final numerical 
results for these corrections. 



2 Detailed description of the model 

For the purpose of calculating the electromagnetic corrections we describe 
low energy n^p scattering by means of individual potentials for each of the 
partial waves which are used in relativised Schrodinger equations (RSEs). We 
therefore develop the model by motivating the construction of the RSEs and 
discussing the components of the potentials that appear in them. 

The point charge Coulomb potential is 

\/P-(r) = a/r, (1) 



where r is the distance between tt"^ and p in the cm. system. When the 
potential (1) is used in the nonrelativistic Schrodinger equation, one obtains 
the nonrelativistic point charge Coulomb amplitude 

fNR = ^^exp{2i((To)jvR-«r7ln(sin2-^)}, (2) 



where 

{cro)NR = aTgT{l + ir]), 



" = ^' <" 

t = -2g2(i_cos^), 



6 being the cm. scattering angle and rric the reduced mass of the n^p system: 

rripfic 



rrir. 



mp + fie' 



with rUp, fie the masses of the proton and charged pion respectively. The 
quantity Qc is the cm. momentum of the 7r~^p system. The Born approximation 
to the point charge Coulomb amplitude is given by the expression (2) without 
the phase factor. 

To make the relativistic generalisation of these nonrelativistic results, we start 
from the one photon exchange (I7-E) no-flip and spin-flip amplitudes fi^E, 
Qi-yE) modifled by the pion and proton form factors F'^ and Ff , Ff respectively, 
and separate them into three parts: 

fl^E = flfE + fl^E + /17E5 (4) 

gifE = gi-yE + 5*17^ + 9vyEy (5) 



where 



and 



fi^E = ^^ , 9T,E = 0, (6) 

fT,'E=^-^iF-F[-l),9ti = 0, (7) 



^ W'-ml-fil , [W - imp - fief] [W^ - imp + fi^f] 



E={ml + qlY/' 



The quantity W is the total energy in the cm. frame and the results in Eqs. 
(6)-(10) are fully relativistic. 

We now note that the relativistic I7-E point charge amplitude /f^^ in Eq.(6) 
is the Born approximation to f^^ multiplied by the factor /c. It is thus the 
Born approximation to the amplitude obtained from the RSE 

{ v' + ql - 2mJ,VP%r) } ^{r) = 0. (11) 



The inclusion of fc in Eq.(ll) gives the unambiguous relativistic generalisation 
of the nonrelativistic potential term 2mcV^'^. The factor /c = 1 when W = fic+ 
rUp and increases as W increases. The RSE (11) is obtained in a more rigorous 
way by following the standard route from the Bethe-Salpeter equation, making 
a three-dimensional reduction, using only the leading part of the I'-fE diagram 
contribution and converting to coordinate space. This is the method used in 
Ref. [10]; the factor fc is just the quantity a in Eqs. (2) and (4) of that paper. 

The full amplitude f^EL obtained from the RSE (11) is 

Irel = ^ exp {2iao - ir]fc In (sin^ -6)}, (12) 



where 



ai = aTgT{l + l + irifc). (13) 



The expression (12) is the point charge Coulomb amplitude to all orders in 
1]. It agrees with the results in Eqs. (2) and (6) in the appropriate limits and 
generalises those results. 

From the RSE (11) one obtains the radial RSEs for individual partial waves: 
( ^ - ^7^ + 1c - 2mJcV^%r) ) ui{r) = 0. (14) 



The wavefunction regular at r = has the asymptotic behaviour as r ^ 00 

sin(gcr - r]fcln{2qcr) - l7i/2 + ai), 

where ai is given by Eq.(13). We shall now proceed to add other potentials to 



V^"^, so that the wavefunction regular at r = behaves as r ^ oo hke 
sm{qcr - r]fc ln(2gcr) - ln/2 + ai + 6i±). 

The notation 6i± takes account of the fact that the extra phase shifts, as well 
as the added potentials, in general depend on the orbital angular momentum 
/ and on the total angular momentum j = ^ ± |. We now take the sum 
over all partial waves to obtain the full no-flip and spin-flip amplitudes /, g 
respectively and use the identity 

^na+5) _ ^ ^ (g2ia _ ^) ^ e^^'^ie^'^ - 1). 

Furthermore, as is customary practice, we remove the phase factor e^*'^° from 
both amplitudes; this does not affect any observable. The result is 



oo 



/ = F' + 5:e2^('^-'^»){(/ + 1)/,+ + lfi^}Pi, (15) 

;=o 

oo 

g = tJ2e''^'''-''^Hfi+-fi-)Pn (16) 



=1 



with 



exp {2i5i±) - 1 
^^^ = 2^^, ' ^^^^ 

/^^= ^^exp{-.r^/,ln(sin2i^)}. (18) 



The potentials other than V^'^ are sufficiently well behaved as r ^ oo for the 
sums over partial waves to be convergent. 

Two photon exchange contributions give a very small correction to the I'-fE 
amplitudes and can be neglected at the level of accuracy required for the 
present calculations. However, even though vacuum polarisation is of higher 
order in a, it needs to be taken into account because it is of extremely long 
range compared with the hadronic interaction. For vacuum polarisation we 
used the standard Uehling potential for point charges given by Durand [11]: 

l^^P(r) = VP%r)I{2mer)2a/37r, (19) 



where me is the electron mass. The explicit integral representation of / can 
be taken from Ref.[ll]. Since V"^ is a potential of electromagnetic origin, 
we assume (as is done in the treatment of pp scattering) that, when V"^ is 



added to V^'^ in the RSEs (14), the factor fc should multiply it as well. The 
tiny effect of the extended charge distributions on V"^ can be neglected. The 
vacuum polarisation amplitude f""^ is given with sufficient accuracy for the 
analysis of current ir^p data by the lowest order expression in Eq.(12.2) of 
Ref.[ll]: 

fp = -^^ {1- cos ey^F (cose), (20) 



where 



i^(cos^) = -^ + X + (l-lx)(l+X)V2ln{ii±|lJ^}, 



4m: 
Ji. — 



t 



The next step is to explicitly separate from the phase shifts 6i±. the very small 
phase shifts crf^^, cr[^' and aj^. The partial-wave projections of the amplitudes 
given in Eqs.(7)-(9) and (20) lead to the results 

+1 
ar' = cxmj^q, J dzt-^Pi{z){FlF^ - 1), (21) 



+1 
c-' = --^^i^ j dzPi(z)F?F^ ± .„... . r^M — ^^ X 
-1 
+1 



-'- ""-^^ dzP,{z)FiF-±^—^ - 

2W J ^ ' ^ 41^(/±l/2 + l/2) 



-1 



22) 



„ _ _ °"lf' jay^l + l,j)(l - ;,)1V'Q<(1 + "V-'). (23) 



a, 







where u = 2ml/q'^ and Qi is the Legendre function of the second kind. The 
expression in Eq.(23) comes from Eq.(8.2) of Ref.[ll]. For / = 0, ...,3 this 
result was checked by direct integration of the RSEs, integrating to 1000 fm 
because of the very long range of V'"^. The expressions (21)-(23) are completely 
sufficient for these very small phase shifts. 

In order to calculate the electromagnetic corrections (which we do for / = 0, 1) 
it is necessary to construct potentials V^f^ and V^o+'i± which reproduce the 



phase shifts ctq^/ and a^^ ^^ with good accuracy up to T^ = 100 MeV. The 
phase shifts themselves are calculated from Eqs.(21) and (22) with the dipole 
form factors 

Ff(t) = (l-t/AJ)-^ 

F|(t) = ^^Ff (t), 
F^t) = {l-t/Alr\ 

where Ap = 805 MeV, A^^ = 1040 MeV. The parameters Ap and A^^ are chosen 
to correspond to the measured charge radii of the proton and charged pion, 
which can be found in Refs. [12,13]. All numerical constants not given explic- 
itly (e.g. the anomalous magnetic moment Hp of the proton) are taken from 
Ref. [14]. The radial RSEs (14) are integrated outwards from the origin, with 
ypc rgpiacgfi Iqj yext Qj, yrei ^^^ ^j^g factor f^ included. 

Even though dipole form factors are used in calculating (Tq^/, we found that 
it is quite sufficient for reproducing these phase shifts and calculating the 
electromagnetic corrections to use the simple potential for gaussian charge 
distributions 

V^''\r) = a/r{eri{r/c) - 1}, (24) 



where 

c^ = 2/3{< r^>p + < r^ >^}. 



Since af^^ in Eq.(21) contains rricfc = (W^ — nip — fil)/2W, it is the potential 
term 2mcfcV'^^^ , with y^^* given by Eq.(24), that reproduces af^* satisfactorily. 

For o"[^' the numerator (W'^ — nip — /i^) of nicfc does not appear explicitly 
in Eq.(22). The potential term can be written as 2mcVf^'' or as 2nicfcV{^\ 
with V{^^ energy independent in each case, and the phase shifts o"q!j_' and a^^f 
satisfactorily fitted. We checked that this makes an insignificant difference to 
the contribution of V{^^ to the electromagnetic corrections. For convenience 
we therefore took the term in the form 2mc/cV^±^ The phase shift cTq'^ is 
proportional to q^ near threshold and can be reproduced very well for T^^ < 100 
MeV by means of a very short range energy independent potential V^^. We 
constructed V^^ and V{1^ by adding together a very short range potential 
and a potential which has the r~^ behaviour for large r characteristic of a 
spin-orbit potential and is modified for small r to take account of the charge 
distributions. These three potentials are plotted in Fig. 1. 
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Fig. 1. The potentials V^^f and V{1^. 
If we now construct the total electromagnetic potentials 



i/em 



ypc j^ yext _^ yrel _^ yvp 



(25) 



and insert them in the radial RSEs of the form given in Eq.(14), with V^'^ 
replaced by VJ^™, the wavefunction regular at r = has the asymptotic be- 
haviour as r ^ oo 



sm 



,r - r]f, ln(2g,r) - /7r/2 + ai + a^' + < + <p). 



Here we have neglected higher order contributions in a to the total electro- 
magnetic phase shifts. The distances at which the asymptotic behaviour is 
reached are determined by V"^ . Since the potentials V[^ have an r~^ tail 
one might think that they require integration of the appropriate RSEs to the 
largest distance. It turns out however that one can stop the integration for 
V{^ at a smaller distance than for V"^. 



We now return to Eqs.(15)-(17) and write 

5i^ = ar + a\i + ar + 6^^. (26) 

The residual phase shifts 5y_, which arise because of the hadronic interaction, 
are usually called nuclear phase shifts, in contrast to the hadronic phase shifts 
to be defined later. We now use the identity 

where 

Ez± = {ai - ao) + ar' + < + a^ (28) 

and 

0"/ — o'o = X! arctan( — ^), (29) 

n=l ^ 

rel 



using Eq.(13). The separation in Eq.(27) is made because the phase shifts af^ 
and o"['^ drop off slowly as / increases and all partial waves must be taken into 
account. It is convenient to treat af"^* in the same way as al^ and a^^ and to 
take the sums indicated in Eqs.(15) and (16) over all values of /. If we work 
consistently to the lowest order in a ( a for ext and re/, ar) for t>p), we can 
replace 



g2i(ai-ao)(^g2i(a-*+a-'+a;^) _ ^^) 



by 



(e^-r' _ 1) + (e^^< _ 1) + (e2-r _ l] 



The sum over all partial waves then yields 

oo 



1=0 



9 = g^E + ^ E(e'^^'Vr+ - e^'^^-fnPl. (31) 



with 



e^^-^ri _ 1 



//+ = • (32) 
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The electromagnetic amplitudes are defined in Eqs.(7)-(9), (18) and (20). The 
electromagnetic phase shifts S^-t are given in Eqs.(28), (29) and (21)-(23). 

We want to emphasise that the expressions given in Eqs.(30) and (31) are 
completely sufficient for the phase shift analysis of the present experimental 
data for T^ < 100 MeV. The lowest energy at which we shall give results is 
Ttt = 10 MeV, for which rjfc = 0.0203. Present experiments do not go to an 
energy as low as this, and the measurements are at angles sufficiently away 
from the forward direction that the lowest order approximation to /^^ given 
in Eq.(20) is adequate. 

For small angles and energies, the leading correction to the lowest order ap- 
proximation to pP is the imaginary part of the quantity a^^^ /qc-, where a\^^ is 
given in Eq. (18.2) of Ref. [11]. The magnitude of the ratio of this correction 
to the lowest order approximation increases as the scattering angle decreases. 
For Tjr = 10, 15 and 20 MeV the ratio reaches 10% at values of 9 very close 
to 3°, 2° and 1° respectively. At T,r = 15 MeV, where 7r"'"p and Ti~p elastic 
scattering are currently being measured, the magnitude of this ratio is 6.6% 
at ^ = 5° and it might just be necessary to take the correction into account 
at such a small angle. 



The final step is to introduce the hadronic interaction via effective hadronic 
potentials VJ^ (in the presence of the electromagnetic interaction) which are 

rem . 
l± ■ 

Vi± = V,T + V.l (33) 



added to V,""" 



Eq.(33) implies that the hadronic potential term has the form 2mc/cV5^. With 
VJ^ chosen to be energy independent, this means that the term has the specific 
energy dependence of the factor fc- It is equally possible to take the hadronic 
potential term as 2mcV/^, with of course a different energy independent poten- 
tial. For the 7r~^p case we checked that the values of the electromagnetic cor- 
rections are practically independent of the energy dependence of the hadronic 
potential term that is chosen. As we shall see in the next paper, this is not the 
case for the two-channel 7r~p situation; the choice of the energy dependence 
of the hadronic potential term has significant consequences for the values of 
the electromagnetic corrections. We shall defer to that paper a full discussion 
of this delicate point of principle; for the present it is sufficient to note that, 
for the 7r~^p corrections, what energy dependence is chosen for the hadronic 
potential term is not important. In the calculation of the corrections given in 
Section 4 we use RSEs with the form 

r/2 1(1-^1) 
' — - ^^ + q! - 2mJ,V,^{r) ) w,±(r) = 0, (34) 



dr'^ r 
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with Vi± given by Eqs.(33) and (25). 

The solutions of the RSEs (34) regular at r = have the asymptotic behaviour 

as r ^ oo 

sin(g,r - r]f, ln(2g,r) - l7i/2 + ai + af^' + < + aj^ + S^^). 



The hadronic phase shifts (5y_ are defined by the asymptotic behaviour as 

r ^ oo 

sin {qcT - In/ 2 + (5^^) 

of the solutions regular at r = of the RSEs in which Vi± in Eq. (34) is replaced 
by V^^. The electromagnetic corrections Ci± are then 

a± = 5-^ - <5f±. (35) 

They arise from the interplay of the potentials VJ^ and Vf^. The quantity Ci± 
can, again neglecting terms of higher order in a, be split into four parts: 

Ci± = Cf^ + Cfl' + Cli + C;"|, (36) 



where each of the parts arises from the interplay of V^ and the correspond- 
ing part of VJj™. We emphasize again that the V^, as well as V^^* and V'^'^\ 
are needed only for the calculation of the electromagnetic corrections. Once 
these corrections are determined, the phase shift analysis, using Eqs.(30)-(32) 
and (35), determines hadronic phase shifts without any further reference to 
potentials. We proceed in the next section to describe the construction of the 
hadronic potentials and the evaluation of the corrections. 



3 Method of evaluation of the corrections 



Numerically we evaluated the electromagnetic corrections Ci± only for the 
three lowest partial waves (0+, 1±). This was done in an iterative procedure 
of which each step consisted of three stages: 

a) from values of the hadronic phase shifts (5/^ we determined hadronic poten- 
tials V^^ which reproduced these 5fj^-^ 

b) using these VJ^ we calculated values for Ci±\ 

c) with these values of Ci± we determined via a phase-shift analysis (PSA) 
new values for the 5y_. 

With these new values of df^ we started again at stage a). The starting phase 



12 



shifts for the whole procedure are of no great importance; in practice we used 
the output from the analysis of Arndt et al. [15]. The iterative procedure was 
continued until the Ci±, and therefore the 5y_, achieved stable values; three 
iteration steps were found to be sufficient. We now describe the three stages 
in more detail. 

a) For each V^ we used a parametric form containing a range parameter 
which was fixed at 1 fm. This form contains three further parameters for each 
partial wave. For a fixed set of these parameters we inserted V^ instead of Vi± 
in Eq.(34) and by numerical integration obtained the solutions u^±{r) regular 
at r = 0. Outside the range of Vl\. the asymptotic form as r ^ oo of the 



Ui^[r} 



ML(r) ~ a\±3i{qj-) + ^Ln^gcr), (37) 



gave the values of (5/^ generated by the VJ^ via the equation 

tan<5L = ^. (3J 



The three parameters in each VJj_ were then varied in order to get the best 
possible fit to the values of (5/j_ up to T^ = 100 MeV. The potentials constructed 
in this way reproduced these values within their experimental errors. In Fig. 
2 we show the potentials V"g'!,_ and V^^ for the final step of the iteration. The 
analytical form and further details can be found in Ref. [16]. 

b) Using the V^^ constructed in a) the electromagnetic corrections were eval- 
uated in the following way. We integrated Eq. (34) with the full potentials Vi± 
and obtained the solutions ui± regular at r = 0, integrating outwards to a 
distance where VJj_ + V^^^ + VJ!j5' + V^^ is negligible compared to V'^'^. In this 
region the asymptotic form of the ui±, 

ui±{r) ~ ai±Fi{r]fc; Qct) + k±Gi{r]fc; Qct), (39) 



where Fi and Gi are the standard point charge Coulomb wavefunctions, yielded 
the nuclear phase shifts S^^ via the equation 

tan (ar* + alt + a7 + 5^^) = ^. (40) 

Eq.(35) then gave the corrections Ci±. for this step of the iteration procedure, 
the 6j[^ being of course the hadronic phase shifts generated by the potentials 
V^^. Note that the use of the asymptotic forms (37) and (39) instead of those 
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Fig. 2. The hadronic potentials Vq^. and V\j^. 

given just before Eq.(35) enables the numerical integration to go to smaller 
distances. 



Though the PSA of the vr+p data is not the topic of this paper we see that 
because of the iterative procedure we have described, the PSA goes hand 
in hand with the evaluation of the electromagnetic corrections. We therefore 
describe briefly how the PSA is done; a full description will be given in a 
separate paper. 

c) The PSA uses Eqs.(30) and (31) with the sums over / taken to / = 3. For 
/ = 2, 3 the small and poorly known nuclear phase shifts (5J^ need to be taken 
from a PSA at higher energies and then extrapolated to the lower energies with 
which we are concerned. Taking the electromagnetic corrections Ci± as calcu- 
lated in stage b), the PSA determined the three hadronic phase shifts (5q_,_, b\j^. 
Parametric forms for these phase shifts were used and the parameters varied 
to obtain the best fit to the data. The parametric forms are simple low energy 
expansions in the pion cm. kinetic energy. Three parameters were used for the 
s-wave and two for the non- resonant p-wave. For the resonant phase shift b\j^^ 
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a parametrised background term was added to a Breit-Wigner resonant term 
with an energy dependent width. The hadronic potentials constructed in a) 
can be considered as an auxihary alternative parametrisation of the hadronic 
phase shifts which is used only for the calculation of the electromagnetic cor- 
rections. 

A very important observation of a conceptual nature needs to be made. We 
have called the quantities Sjj_ hadronic phase shifts, but they are obtained 
from a phase shift analysis which contains the physical masses of 7r+ and p. 
It is therefore clear that the (5[^ cannot be considered to be strictly hadronic 
quantities. Such quantities relate to a situation in which the electromagnetic 
interaction is switched off {a = 0) and therefore need to be obtained using the 
hadronic masses of tt"*" and p. These hadronic masses are not the same as their 
physical masses (indeed it is universally accepted that the hadronic mass of 



TT"'" is very close to the physical mass of 7r°). In the same sense we have noted 
that the quantities V/^. are effective hadronic potentials in the presence of the 
electromagnetic interaction. They would be different in the complete absence 
of all electromagnetic interactions. Therefore in our present work, as in all 
previous work known to us, we avoid speculation about the strictly hadronic 
situation and give quantities that we call hadronic phase shifts and electro- 
magnetic corrections which, though they have a precise definition within the 
framework of our potential model, are not the full corrections that would give 
truly hadronic phase shifts by subtraction from the nuclear phase shifts. Any 
attempt to completely purge the experimental data of all electromagnetic ef- 
fects would be very speculative. Our aim is the more modest one of taking 
account of those electromagnetic effects that may be calculated with reason- 
able confidence. 



4 Numerical results for the corrections 



The final results for the three electromagnetic corrections Co+, Ci± are given 
in Table 1 from 10 to 100 MeV pion lab kinetic energy. Estimates of the 
uncertainties in the values of these corrections are also given in the table. In 
making these estimates it is necessary to distinguish between a) uncertainties 
coming from applying our model and b) uncertainties coming from the choice 
of the model itself. 

One source of a) is the uncertainty in the hadronic phase shifts, which comes 
from the PSA of the experimental data (experimental errors in the data and 
uncertainty in the d- and /-wave phase shifts used as input). The uncertain- 
ties in the hadronic phase shifts will be given in full in the separate paper on 
the PSA that was mentioned earlier. The resulting uncertainties in the cor- 
rections were thoroughly studied by keeping track of the corrections obtained 
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Table 1 

Values in degrees of the electromagnetic corrections Co+, Ci_ and Ci+ as functions 

of the pion lab kinetic energy T-,^ (in MeV) . 



T. 


Co+ 


Ci+ 


Ci- 


10 


0.081± 0.003 


-0.023± 0.000 


0.005± 0.000 


15 


0.083± 0.003 


-0.034± 0.000 


0.007± 0.000 


20 


0.085± 0.004 


-0.047± 0.001 


0.009± 0.000 


25 


0.088± 0.005 


-0.060± 0.001 


O.OlOi 0.000 


30 


0.090± 0.005 


-0.074± 0.002 


0.012± 0.000 


35 


0.093± 0.006 


-0.090± 0.002 


0.013± 0.001 


40 


0.096± 0.006 


-0.106± 0.003 


0.014± 0.001 


45 


0.099± 0.007 


-0.125± 0.003 


0.015± 0.001 


50 


O.lOlib 0.007 


-0.145± 0.004 


0.016± 0.001 


55 


0.104± 0.008 


-0.168± 0.005 


0.017± 0.002 


60 


0.107± 0.009 


-0.194± 0.006 


0.018± 0.002 


65 


0.109± 0.010 


-0.223± 0.008 


0.018± 0.002 


70 


O.lllib 0.010 


-0.255± 0.011 


0.019± 0.002 


75 


0.114± 0.011 


-0.291± 0.016 


0.019± 0.003 


80 


0.115± 0.012 


-0.332± 0.023 


0.019± 0.003 


85 


0.118± 0.014 


-0.378± 0.032 


0.019± 0.004 


90 


0.119± 0.015 


-0.429± 0.045 


0.019± 0.004 


95 


0.121± 0.016 


-0.485± 0.062 


0.019± 0.004 


100 


0.123± 0.018 


-0.547± 0.088 


0.020± 0.005 



using a variety of input values for the hadronic phase shifts (in particular for 
the successive iteration steps described in Section 3). These uncertainties are 
substantially smaller than those coming from the use of the particular form of 
the parametrised hadronic potentials, which are also of type a). For the final 
numerical calculations we fixed the range parameter in the potentials at 1 fm. 
Its value can be varied from 0.8 fm to 1.2 fm without any significant change 
in the quality of the fits to the hadronic phase shifts. In this way we obtained 
an estimate of the uncertainties coming from the fact that some of the fine 
details of the potentials may have been missed. 

The estimated uncertainties given in Table 1 were obtained by combining the 
two sources of type a) described in the previous paragraph. For Co+ and Ci_ 
these uncertainties are much smaller than the errors on the hadronic phase 
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Table 2 

Values in degrees of the various contributions to the electromagnetic corrections 

Co+ as functions of the pion lab kinetic energy T^ (in MeV) . 



T^ 




r~iext 
'-'0+ 




'-'0+ 


10 


0.093 


-0.008 


-0.002 


0.000 


20 


0.102 


-0.013 


-0.004 


0.000 


30 


0.113 


-0.018 


-0.005 


0.000 


40 


0.123 


-0.022 


-0.005 


0.000 


50 


0.134 


-0.027 


-0.006 


0.000 


60 


0.145 


-0.031 


-0.007 


0.000 


70 


0.155 


-0.036 


-0.008 


0.000 


80 


0.165 


-0.041 


-0.009 


0.000 


90 


0.175 


-0.046 


-0.010 


0.000 


100 


0.183 


-0.051 


-0.010 


0.000 



shifts that arise from the PSA itself. For example, at 100 MeV the errors 
on 5o+ and 5i_ are 0.15° and 0.09°, compared with the uncertainties in the 
corrections of 0.018° and 0.005° respectively. The situation is different for Ci+, 
where the error on (5i+ at 100 MeV is 0.05°, while the correction is —0.547° 
and its uncertainty is 0.088°. The correction itself is very large (5i+ itself is 
21.27° at this energy) and the uncertainty in the correction is somewhat larger 
than the error arising from the PSA. 

The corrections in Table 1 are intended for use in future PSAs of experimental 
data. Such PSAs need to use Eqs (30)- (32), (28) and (35) as well as higher 
order corrections to /^^ and o^^ for experiments at very low energies and 
angles, as discussed in Section 2. For the output phase shift (5i+ the uncertainty 
in the correction Ci+ needs to be taken into account as well as the statistical 
error arising from the PSA itself. 

Uncertainties of type b) arise from the choice of the RSEs (34), with the 
hadronic potential term having the specific energy dependence of the factor 
/c. We have already remarked that the choice of the energy dependence of this 
term has no significant effect on the corrections. The RSEs themselves are the 
only relativistic equations that are well suited to a two-body problem. It is 
not meaningful to give an uncertainty arising from this source; all one can do 
is to compare the results obtained using our present model with those from 
alternative models, as we do in a moment. 

In Tables 2 and 3 we give (at a smaller set of energies) the various contributions 
to the corrections (Eq.(ll)). These components are additive to a very good 
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Table 3 

Values in degrees of the various contributions to the electromagnetic corrections 

Ci_ and Ci-i- as functions of the pion lab kinetic energy T^ (in MeV). 



T. 


cll 


r^ext 


^1+ 


cz 


10 


-0.023 


0.000 


0.000 


-0.000 


20 


-0.047 


0.001 


-0.000 


-0.000 


30 


-0.074 


0.001 


-0.001 


-0.001 


40 


-0.104 


0.002 


-0.004 


-0.001 


50 


-0.141 


0.004 


-0.008 


-0.001 


60 


-0.184 


0.007 


-0.016 


-0.001 


70 


-0.236 


0.011 


-0.029 


-0.002 


80 


-0.300 


0.017 


-0.048 


-0.002 


90 


-0.376 


0.025 


-0.077 


-0.003 


100 


-0.464 


0.036 


-0.119 


-0.003 


T. 


clI 


(jexj 


CI"} 


cz 


10 


0.004 


-0.000 


0.000 


0.000 


20 


0.008 


-0.000 


0.000 


0.000 


30 


0.011 


-0.000 


0.000 


0.000 


40 


0.014 


-0.000 


0.000 


0.000 


50 


0.016 


-0.000 


0.000 


0.000 


60 


0.018 


-0.000 


-0.000 


0.000 


70 


0.019 


-0.000 


-0.001 


0.000 


80 


0.021 


-0.001 


-0.001 


0.000 


90 


0.022 


-0.001 


-0.002 


0.000 


100 


0.024 


-0.001 


-0.004 


0.000 



approximation and tiny differences between the sums of the numbers in these 
tables and the complete corrections in Table 1 are due to higher order effects. 
The corrections are clearly dominated by the component C^'^ arising from the 
interplay between the hadronic and point charge Coulomb potentials. The 
only other significant components are Cg^ and C[^' near T^ = 100 MeV. 



Our results are compared with those of the dispersion theory approach used 
by N0RDITA[3] in Fig. 3. They have not included the vacuum polarisation 
contribution but we see from Tables 2 and 3 that it is negligible for each of 
the partial waves at the energies we consider. The results of Zimmermann [8], 
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Fig. 3. Values in degrees of the electromagnetic corrections Co+, Ci+ and Ci- from 
our present calculation (solid curves), from NORDITA [3] (circles) and Zimmermann 
[8] (triangles). 



who uses a sUghtly different version of the potential model, are also shown in 
Fig. 3. We have added our values of the relativistic corrections to the results of 
Zimmermann. We have indicated in Fig. 3 the uncertainties in our corrections 
at 100 MeV, as given in Table 1. No errors are given in Refs. [3] and [8] for the 
corrections presented there. For Co+ the results are in very good agreement, 
the differences being considerably smaller than the error on 6q_^_ arising from 
the PSA. For Ci+ our results agree quite well with those of Zimmermann, 
which indicates the stability of the results obtained with the potential model. 
However, there is a systematic difference from the results of NORDITA, the 
discrepancy increasing with energy. At 100 MeV the difference is about 2.5 
times our estimated uncertainty in Ci+ and over 4 times the error in 5i+ 
arising from the PSA. For Ci_ our results agree with those of NORDITA at 
low energies, but at 100 MeV the discrepancy is roughly twice the error in 
(5{'_ . However, since S'l_ is very small, the corrections are of minor importance 
for the results of the PSA. 
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As we discussed in Section 1, the calculation of electromagnetic corrections us- 
ing dispersion relations omits what could be important medium range effects 
due to t-and -u-channel exhanges, which the potential model includes quite re- 
liably. The differences of our results from those of NORDITA are probably due 
to such medium range effects. In particular the stronger energy dependence of 
the NORDITA results for Ci+ shown in Fig. 3 is likely to be a medium range 
effect, perhaps due to their omission of t-channel tttt (T = 0, J = 0) exchange. 
We therefore claim a higher degree of reliability for our present calculation 
of the electromagnetic corrections, compared with that of NORDITA. At the 
same time it is pleasing that the corrections in the n^p case are very reli- 
ably known and are to a large extent independent of the model used for their 
calculation. 
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